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We consider phase separation in nonequilibrium Bose gas with an attractive interaction between 
the particles. Using numerical integrations on a lattice, we show that the system evolves into a state 
that contains drops of Bose-Einstein condensate suspended in uncondensed gas. When the initial 
gas is sufficiently rarefied, the rate of formation of this quantum dew scales with the initial density 
as expected for a process governed by two-particle collisions. 
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Theory of interacting Bose gases has been an impor- 
tant part of quantum statistical mechanics ever since Bo- 
goliubov's seminal work |l|. In particular, nonequilibrium 
Bose gases are interesting from at least two points of 
view. First, such gases can now be produced in labora- 
tory via modern cooling techniques. A dramatic demon- 
stration of the resulting nonlinear dynamics is Bose- 
Einstein condensation (BEC) observed in alkali vapors 
P-H. Second, nonequilibrium gases of elementary parti- 
cles frequently arise in cosmological scenarios and could 
have played an important role in the evolution of the 
universe. 

One way how nonequilibrium Bose gases arise in cos- 
mology is via decay of a coherently oscillating field. This 
mechanism could be important, for instance, at the end 
of an inflationary stage, i.e. during the reheating af- 
ter inflation. Indeed, it has been found that in some 
inflationary models, the oscillating inflaton field decays 
rapidly and completely into a gas that contains both 
the inflaton quanta and other types of Bose particles 
[g|. These gases have very large occupation numbers in 
low-momentum modes and almost no occupation in high- 
momentum modes; they are highly nonthermal. The pos- 
sibility of existence of such gases is by no means limited 
to the postinflationary era. In particular, there are in- 
dications that nonbaryonic cold dark matter constitutes 
a significant fraction of the matter in the universe at 
present. At the epoch of galaxy formation, gravitational 
instability develops on a variety of scales, which may lead 
to formation of small-scale dark matter clumps in galaxy 
halos. Dark matter particles trapped in a gravitational 
well are out of thermal equilibrium and are nonrelativis- 
tic. Typically interparticle interactions are very small, 
but, if the particles are bosons, the relaxation time can 
in certain cases be comparable to the age of the universe 
Q . This opens a possibility of Bose-Einstein condensa- 
tion and formation of Bose-stars m-M- One proposed 
precursor of those is axion miniclusters |10[| , but modern 
particle models contain a variety of other fields of poten- 
tial interest in this respect: majoron, dilaton, moduli, to 



name a few. 

Thus, it is important to investigate the evolution of 
nonequilibrium Bose gases under various conditions. If 
the interaction between the particles is repulsive, and the 
energy density is sufficiently low, a Bose-Einstein con- 
densate will form. The process of Bose-Einstein conden- 
sation in this case has been studied theoretically in a 
number of papers ||] , |ll|-|l3| . 

The question we want to address in this paper is what 
happens to a nonequilibrium Bose gas if the interaction 
between its particles is attractive, at least within a cer- 
tain range of interparticle distances. There is hardly any 
doubt that an attractive interaction will lead to clump- 
ing and phase separation, and statements to that effect 
have appeared in recent literature [H. However, it has 
remained unclear whether the clumps will be in the nor- 
mal or the superfluid state. In addition, kinetics of the 
clumping needs to be elucidated. Our main result, ob- 
tained via numerical integrations, is that the clumps are 
drops of Bose-Einstein condensate, i.e. each of them is 
characterized by a macroscopic order parameter. These 
drops remained suspended in uncondensed gas for as long 
as we could follow the evolution, although they did grow 
somewhat at the expense of the gas. Because Bose con- 
densation in the drops is attributable to the quantum 
statistics of the particles, we call such drops quantum 
dew. 

The coherent, macroscopically ordered nature of quan- 
tum dew may be important in cosmological (as well as 
laboratory) applications. Suppose for example that the 
particles it is made of can decay into some other parti- 
cles. The macroscopically populated mode of a coherent 
clump may work as a laser P,15[; as a result, quantum 
dew may decay much faster than an incoherent clump 
would. 

The purpose of the present paper is to prove the co- 
herent nature of the clumps and to study the kinetics 
of appearance and growth of dew drops. For this pur- 
pose, we have chosen the simplest model with an attrac- 
tive interaction and, nevertheless, a stable ground state. 



From the point of view of cosmological applications, per- 
haps the most important effects left out of this simple 
model are the expansion of the universe and the gravita- 
tional attraction. The expansion dilutes gas available for 
clumping and thus slows the clumping down. The gravi- 
tational attractions works in the opposite direction. The 
net effect of these opposing tendencies can in principle be 
found via numerical integrations, and we plan to return 
to this important question in future. 

The model contains a nonrelativistic complex Bose 
field tp with the following equation of motion 



at 
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The field is normalized so that the attractive cubic term 
on the right-hand side has the coefficient of unity. The 
corresponding coupling 174 then appears in the commuta- 
tion relation: 



[ak'«k'] = Iff4|'^kk' ' 



(2) 



the annihilation operators aj^ being defined via "0(1") — 
y-i/2 ^ ^j^ exp(ikr) in a finite volume V . The quintic 
term in p) is repulsive, and it becomes important when 
'0t'0 approaches g^ . (The coupling 54 appearing in (g) 
is related to A of the relativistic A0'*/4 potential via 54 = 
3A/27TI, and to the scattering length a of a nonrelativistic 
Bose gas via 54 = Stto; here ?i = 1. The physical density 
is ■>pm\g4:\.) 

Our integrations are set up as follows. In the initial 
state the occupation numbers nj^ 
sian distribution over momenta 



'^If'^k have a Gaus- 



^exp(-fcVfc2) 



(3) 



The population of the homogeneous mode is no — A, 
which is not considerably larger than population of other 
modes with small k. In this sense, there is no macroscopic 
condensate in the initial state. 

Now, we assume that (74 in (0) is small compared to A 
in (J^). Then, we can neglect the commutator of a and 
a^ compared to the typical magnitude of a itself, cf. Bo- 
goliubov y . As a result, the problem becomes classical 
and can be integrated on a lattice. This classical approx- 
imation has been used to study nonlinear dynamics of 
relativistic Bose fields at large occupation numbers fll] 
and the process of Bose- Einstein condensation |1^. In 
the present work, we use it to demonstrate the forma- 
tion of quantum dew in the model p) . This use involves 
no contradiction of terms: quantum dew is an effect of 
quantum statistics when one thinks in terms of individual 
particles, but it comes out as an effect of classical evolu- 
tion in the collective, field-theoretical description. Notice 
that eq. (|^) determines only absolute values of aj^; their 
phases are chosen as uncorrelated random numbers. 

We choose the parameters of the model in such a way 
that {ip^ip) <gC gQ^; the angular brackets denote averag- 
ing over the lattice. This means that in the initial state 



the attractive interaction is much more important than 
the repulsive one. In this case, we expect that, in appro- 
priate dimensionless units, the time scale tc of the initial 
collapse of the gas into clumps depends only on the single 
remaining parameter of nonlinearity 



e = 



2me 



(4) 



where e is the average kinetic energy per particle in the 
initial state; ^ is of order of the ratio of the initial po- 
tential energy of attraction to the initial kinetic energy. 
A similar parameter for an atomic gas in a trap will be 
introduced below. We can write 



tc'=^FiO 



(5) 



where F is some function obeying the condition F(0) = 
and F{1) ~ 1. The form of F(^) at small ^ is established 
below. For the initial distribution (p|), the initial param- 
eter of nonlinearity is ^ = koA/12Tr'/^. We choose units 
of time so that 2ni — 1. We also use units of length in 
which /cq = 27r, i.e. measure lengths in units of the parti- 
cles' typical initial de Broglie wavelength. Except where 
stated otherwise, we consider the case of moderate non- 
linearity A ^ 5, which corresponds to ^ = 0.47. We use 
g-i = 3600. 

The results below are from integrations on a 64'^ cubic 
lattice with side L = 2.25 (in the above length units) and 
periodic boundary conditions. The state of the system 
was updated via a second-order in time algorithm based 
on the Crank-Nicholson method for the diffusion equa- 
tion. The algorithm conserves the number of particles 
exactly. Energy non-conservation was below 2% for the 
entire integration time. 

Fig. 1 shows two snapshots of the field, at times i = 0.1 
and t — 20. Dots have been placed on all lattice sites at 
which 1-01 > 30 (the mean-square value of |0| is 5.3). 
Drops of dew are clearly seen. 




FIG. 1. Drops of dew at different moments of time. 

A movie of the evolution of this picture from i = 0.1 to 
t = 20 shows that the drops of dew move around, gradu- 
ally slowing down, and occasionally coalesce. The over- 
all growth of the number of sites with \tp\ > 30 continues 
even at i = 20, the latest time in our computation, but at 
that time it is already quite slow. If we define that grid 
points with \^\ < 8 belong to the gas, and correspond- 
ingly grid points with j^l > 8 belong to the dew (j^l — 8 



is approximately the boundary between the gas and the 
dew at t = 20, see Fig. y below), we find that around 15% 
of all particles are in the gas, and around 85% had con- 
densed in the dew by the time t — 20. These fractions, 
however, may be altered when gravity is included. 
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FIG. 2. Condensation process. 
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Fig. shows initial stages of the condensation process: 
we plot the fraction of particles that are in the dew, as a 
function of time. We observe two distinct stages: a rapid 
collapse followed by a slower chaotic evolution. Because 
f ^ 1 at ^ = 5, we estimate the time of collapse tc from 
(a) as tc ~ ^q"^ ■ For ^o = 27r this gives tc ^ 0.025, in 
good agreement with the data of Fig. g. In the regime 
of weak nonlinearity, ^ <C 1, we expect the collapse to 
be due to two-particle collisions, in which case F(^) ex 
^2. (With this form of F(^), the estimate (|) for the 
time of the collapse coincides with the estimate of the 
condensation time of Ref. [pj , which can also be obtained 
from a solution to the Boltzmann equation |E2|). The 
time of the collapse then has to scale as A^^^hen we 
decrease A and keep all other parameters fixed. Results 
of integrations with different values of ^ < 1 confirm this, 
see Fig. |[ 

For an atomic gas confined in a trap, at some temper- 
ature T, one can introduce initial parameter of nonlin- 
earity ^p: 



Ct 



47r?i \a\n 
rriT 



(6) 



where n is a typical gas density, and a is the scatter- 
ing length, which in the present case is negative. In (g), 
fcs = 1, but we have restored fi. Let us use for estimates 
n = (mro/3.31?i^)3/2 ^nd T = To, where To is the tem- 
perature of BEC of an ideal monoatomic gas in a given 
trap and with a given number of particles. Bradley et 
al. quote Tq = 300 nK and a = — 27.3ao for their ex- 
periment with trapped ^Li (ao is the Bohr radius); using 
these values we obtain ^t — 0.006. Estimating the rate 
of the collapse as ?iij^ ~ ^oCt, we find ic ~ 1 s. We thus 
expect that quantum dew can be observed in laboratory 
in traps of a sufficiently large size. 



The onset of the slower chaotic evolution indicates that 
a chemical quasiequilibrium between the dew and the gas 
has been reached, i.e. the processes of evaporation of par- 
ticles from the existing dew drops and condensation back 
onto them are approximately (but not exactly) balanced. 
This interpretation is supported by the following test. 
The probability distribution of the absolute value of the 
field over lattice sites shows two distinct peaks: one at 
large 1-01, corresponding to the dew drops, and another at 
small IV'I, corresponding to the gas of particles, see Fig. 
3. If at some instant we remove the gas, i.e. set i/" = 
at all sites where we had \ip\ < 10, and then continue the 
evolution, the gas reappears, while the number of sites 
occupied by the dew decreases down to another slowly 
evolving value. Apparently, the dew partially evaporates, 
so as to restore the chemical quasiequilibrium. 







FIG. 3. Probability distribution function of the field mag- 
nitude at t = 20 (the solid line). Initial probability distribu- 
tion is shown by the dotted line. 
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FIG. 4. Field distribution in a spatial slice at t = 0.2. 



Finally, Fig. 4 illustrates the coherent nature of the 
dew. It shows the field tp at t = 0.2 at a section of our 
integration cube parallel to the x-y plane. For visual 
clarity only even-even numbered sites are included. The 
length of an arrow represents [ipl, and the angle clock- 
wise from 12 noon represents aigip. We see that the 
sites occupied by the dew (i.e. having large |V'|) are in 
drops, and each such drop is coherent — the arrows point 
approximately in the same direction. The direction of 
arrows in each drop rotates with time, just as in the ho- 
mogeneous case 0], but these directions are different for 
different drops. Similar slices at later times show that 
the clumping becomes more pronounced, the dew is still 
coherent, while the remaining gas (occupying sites with 
small \tlA) is incoherent. 

Eq. m) has stable nontopological solitons of the form 



■i(j{r, t) — x(r) exp(iu;i) 



(7) 



(nonrelativistic analogs of Q-balls ||l^). As we continue 
to truncate the gas, i.e. to remove particles from sites 
with progressively smaller ji/;!, and to evolve the system 
between these truncations, we expect to eventually reach 
a state in which solitons float in vacuum (or gas of a 
very small density) . Changes in the probability distribu- 
tion function (p.d.f^ of \'ip\ resulting from this procedure 
are shown in Fig. H. We interpret the limiting form to 
which the p.d.f. converges in the middle range of li/)! as 
corresponding to the wall profile of nontopological soli- 
tons. Computation of ip at the center of a soliton in the 
thin-wall approximation gives \ip\c = (3/4.96)^^^ ~ 52, 
in good agreement with the position of the peak in the 
p.d.f. at large \^\. Like nontopological sohtons produced 
in a decay of an unstable homogeneous condensate [Q , 
quantum dew may work as cold dark matter. 
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FIG. 5. P.d.f of the field's magnitude after truncating the 
gas and letting it reequilibrate. Solid line: no truncation; dot- 
ted line: after truncation at ?/> = 10; dashed line: after further 
truncation at ^ = 3; long dashed: after further truncation at 

7/) = 1. 

To summarize, our main results are: (i) a numerical 
proof that the clumps of matter formed in a nonequi- 
librium gas with an attractive interaction are coherent 



drops of Bose-Einstein condensate; (ii) evidence that the 
rapid collapse of particles into drops of this quantum dew 
is followed by a slower evolution, during which the dew is 
in approximate chemical equilibrium with the surround- 
ing gas; (iii) evidence that at weak nonlinearity the rate 
of the initial collapse is consistent with being determined 
by two-particle collisions. 
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